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Abstract 

We propose a method for extracting the shear viscosity over entropy density ratio (r]/s) of the 
quark-gluon plasma from experimental data. We argue that uncertainty due to poor knowledge 
of the earliest stages of a heavy-ion collision is smallest for ultra-central events. The most pre- 
cise value of rj/s can thus be obtained from a global fit to pj -integrated Fourier harmonics of 
azimuthal correlations. We further outline a method for quantifying the overall uncertainty in 
the extracted value. Only after a comprehensive and systematic accounting of all sources of 
uncertainty can a reliable measurement be claimed. In these proceedings we report preliminary 
results; full and final results will be presented in a separate publication. 



1. Introduction 

Transport coefficients such as shear viscosity represent fundamental properties of a system 
and contain direct information about the underlying physics. For QCD matter in the temperature 
range probed by an ultra-relativistic heavy-ion collision, such quantities are exceptionally diffi- 
cult to calculate from first principles. On the other hand, observables in these experiments can be 
quite sensitive to transport coefficients, potentially allowing for a measurement to be extracted 
from data [Q]]. Thus, such a precision extraction is a topic of significant interest. 

In the past, a precision extraction has been hampered by significant theoretical uncertainty, 
especially that arising from poor knowledge of particle production and thermalization in the 
earliest stages. A number of viscous hydrodynamic calculations have been performed, which 
have been quite successful at reproducing flow measurements J^Sl- Typically, a small value of 
viscosity is required in order to fit data — not too much more than a few times the conjectured 
lower bound r\js > 1 I An. However, the precise best-fit value depends on various aspects of the 
calculations that are not fully constrained, and this introduces uncertainty in any extraction from 
data. Although these uncertainties have been the target of significant investigation, it has so far 
been done in a piecemeal fashion [4]. No one has yet attempted to compile this information 
in a systematic study of all sources of uncertainty in order to place a robust and reliable upper 
bound on rj/s. In addition, the conjectured bound is now known not to be universal, so it is of 
great interest to determine if one can place a lower bound on the quantity. Currently, there is no 
convincing evidence that measured data can not be reproduced with ideal hydrodynamics, and in 
fact most flow data can be reproduced quite well in such calculations fl5[]. In these proceedings 
we present a proposal for how to extract rj/ s with the smallest uncertainty currently possible, as 
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Figure 1 : (Color online) RMS eccentricities versus centrality for several models of initial conditions in 1 % centrality bins. 
Centrality was selected according to multiplicity or entropy, as appropriate for each model, ancRhe transverse density of 
the same quantity was used to calculate eccentricities. Inset plots magnify the most central 3% of collisions. The models 
shown are: PHOBOS MC-Glauber (<| [Glauber (points)], MC-KLN Ql v3.43 and its Glauber implementation [Glauber 
(discs)], mckt-vl.25 with negative binomial fluctuations OJ (MCrcBK (NB)), DIPSY @|, and UrQMD (hJ. 



well as preliminary results on such an extraction including a quantitative error bar obtained by 
systematic study of all known sources of uncertainty. 

2. Ultra-central collisions 

In a standard picture of a heavy-ion collision, the asymmetric region of overlap in a collision 
between two nuclei results in a anisotropy in the transverse density of the system at early times. 
Over time, this spatial anisotropy is transformed into a momentum anisotropy by a collective 
response that is well described by (viscous) hydrodynamics. In fact, the elliptic flow coefficient 
V2 is (to a good approximation) directly proportional to this initial eccentricity in a given colli- 
sion event: \>2 cc g 2 = |(r 2 e 2 "^)|/(r 2 ). The proportionality constant depends sensitively on the 
viscosity of the medium, However, poor knowledge of the earliest stages of the collision result 
in significant uncertainty in so See, for example, the top left panel of Fig.Q] which shows root- 
mean-square values of from a number of commonly-used models for initial condition^. The 
resulting uncertainty in v 2 has been the largest source of uncertainty when determining r\js from 
data J2]. Though there has been much recent study concerning the physics of the early stages of 
a collision fll^l, it is not yet possible to make a reliable statement about what range of values of 
initial eccentricity are within the realm of possibility. 

It is now known that quantum fluctuations are also a significant source of anisotropy in the 
initial state, and therefore anisotropic flow. This is especially apparent from the presence of large 



Most v„ measurements are sensitive to the rms value of the event-by-event distribution of v„ 
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Figure 2: (Color online) Left: Example plot of chi-squared per degree of freedom for a one-parameter simultaneous fit of 
integrated v„ forn = 2-6 at 0-1% centrality, as a function of the ratio of shear viscosity to entropy density r]/s. Right: Fit 
values corresponding to the minimum chi-squared curves in the left plot. The green bars represent measured values from 
the ATLAS collaboration (l3|l; the vertical thickness represents the combined statistical and systematic error. Models 
with large fluctuations and a correspondingly large best-fit r//s tend to have a steeper slope than models with smaller 
eccentricities. The range of best-fit viscosity values for the curves in the left plot gives an indication of the uncertainty in 
an extraction oir]/s due to the uncertain initial anisotropy in the initial state (~ ± 0.05). 



odd harmonics, such as triangular flow V3, which is generated from a spatial "triangularity", de- 
fined similarly to the eccentricity V3 oc g 3 This would be zero at midrapidity in the absence of 
event-by-event fluctuations. The triangularity that one sees in a calculation is largely determined 
by the strength of these fluctuations. In contrast, the differences in particle production mecha- 
nisms that lead to very different £2 in non-central collisions has very little effect on £3. Compare, 
for example, the triangularities of the MC-KLN model versus Glauber (discs) in Fig. [1] which 
contain only fluctuations from random position of nucleons in the colliding nuclei, implemented 
in exactly the same way in the MC-KLN code. The triangularity £3 is very similar, despite the 
fact that the two pictures of particle production result in a very different eccentricity £2. 

By comparing the results from other models, it is clear that there still remains some uncer- 
tainty due to different way fluctuations are implemented. Despite the seemingly modest reduction 
in uncertainty (as compared to the uncertainty in £2), one can reasonably argue that the results 
in the second panel of Fig. Q] represent the full range of results that can be expected, in marked 
contrast to the results in the first panel. The models with the lowest triangularity are known to 
be missing important sources of fluctuations, while the large fluctuations in the models with the 
largest £3 cannot be reconciled with constraints from v\ measurements [11]. 

One might thus propose to extract rj/s from measurements of V3 to minimize uncertainty from 
the initial conditions. However, we can make use of more experimental constraints by noting 
that in the most central collisions, all harmonics are generated entirely from fluctuations, and the 
above discussion still applies. One further convenient factor is that, for ultra-central collisions, 
the approximate proportionality v„ oc s„ holds for all n Jl2ll . This allows for a comprehensive 
parameter study with reasonable computing resources. 

Thus, our proposal is to extract rj/s from a simultaneous fit to pj -integrated v„ measurements 
from ultra-central collisions at the highest-energy heavy-ion collisions. These data are least 
sensitive to uncertain aspects of theory, and therefore allow for the the most precise extraction of 
77/ s with the smallest possible uncertainty. Such data are available from the ATLAS collaboration 
for V2-V6 in the 1% of most central collisions 1 1311 . 



3 



3. Extracting t//s 



In a theoretical calculation, the viscosity of the medium can be taken to be a free parameter. 
With all other parameters fixed, one can calculate v„ as a function of n/ s, and perform a least- 
squares fit to find the preferred value. If a particular parameter is uncertain, one can vary that 
parameter within the allowed range and redo the fit. The shift of best-fit 77/s indicates the size 
of the resulting error bar in 77/5 due to this particular source of uncertainty. If one then varies all 
possible aspects of the theory in all possible combinations, the smallest and the largest extracted 
values thus represent a robust, if conservative, upper and lower error bar for the quantity. 

This is the idea behind the preliminary results shown at this conference. Details will be 
given in a forthcoming publication, but an example fit for a calculation using a single set of 
hydrodynamic parameters but different initial conditions is shown in Fig. [2] It is notable that 
hydrodynamics can, in fact, simultaneously give a good description of all v„. However, it is 
also interesting to note the tendency for tension in the relative values of v% and V3 for all initial 
conditions tested. 

We reported a preliminary result of 0.07 < 77/s < 0.43, but a few comments are in order. First, 
the calculations used a constant value of rj/s during the hydrodynamic evolution. Therefore, the 
result is an average value that corresponds to the viscosity at the typical time when flow is built 
up. The lower bound indicates that purely ideal hydrodynamic calculations can not describe 
these data. This is a significant result. However, it is technically possible that the viscosity could 
dip below our lower bound at some point during the evolution, so this result does not represent 
a lower bound for the viscosity of QCD over its entire temperature range. Second, since it has 
never before been done, we want to be confident at every stage that we have a robust and reliable 
error bar. The preliminary result, therefore, is fairly conservative, and the error bar will likely be 
smaller in the final result after further study. 



4. Conclusions 

We argue that at the current time, the most precise value for the shear viscosity to entropy 
density ratio tj/s of the quark-gluon plasma can be obtained by a fit to /^-integrated v„ mea- 
surements in ultra-central collisions. We further report preliminary results such an extraction, 
including an error bar from a systematic study of all known sources of uncertainty. 

We thank the PHOBOS Collaboration, Adrian Dumitru, Yasushi Nara, and Christoffer Flens- 
burg for providing their code to the public, which allowed us to calculate the eccentricities in 
various Montre-Carlo models, and Hannah Petersen for calculating them for us using UrQMD. 
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